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The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase in the vortex lattice state is quantitatively studied 
using the self-consistent Eilenberger theory in three-dimensional (3D) space. We estimate free energy 
to determine the FFLO phase diagram in the H-T plane and stable FFLO wave number in the isotropic 
system with the 3D Fermi sphere and s-wave pairing. To facilitate the experimental identification of the 
FFLO state, we investigate the field evolution of NMR spectra and flux line lattice form factors obtained 
in neutron scattering in the FFLO vortex states. Possible applications of our results to experimental data 
on CeCoIn 5 are mentioned. 
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There has been much attention focused on discovering 
and hunting exotic superconducting states. Among them, 
in the singlet pairing category, the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state has been one of the most elu- 
sive examples since its theoretical prediction in 1964. 1,2) In 
the FFLO states, the superconducting order parameter ex- 
hibits a spatial modulation. 3) The population imbalance is 
brought about either by the application of an external field 
in the charged electron case through the Zeeman effect or by 
the preparation of up and down species in cold neutral atom 
gases. 4_6) Under an applied field, the FFLO state is expected 
to be the one most likely to emerge in the low- temperature 
(7) high-field (H) region. So far there is no direct evidence to 
prove the existence of the FFLO state in either the charged 
or neutral system. Thus the FFLO state still remains elusive. 
In the charged system, the orbital depairing effect due to elec- 
tron diamagnetic motion in a magnetic field cannot be ignored 
since it may affect the stability of the FFLO state. 

On the other hand, there has been no microscopic cal- 
culation of the field evolution of the FFLO state that fully 
takes account of the orbital depairing effect, namely, the vor- 
tex effect beyond the Ginzburg-Landau (GL) framework valid 
near the upper critical field (H C 2)^ In particular, the Larkin- 
Ovchinnikov (LO) state with a periodically modulated am- 
plitude of the order parameter is highly difficult to describe 
owing to the solitonic spatial variation with infinitely many 
higher harmonics in general. 8) In CeCoIns, it is suggested that 
the LO state is realized rather than the Fulde-Ferrell (FF) state 
in which only the phase is modulated in the order parame- 
ter. 3) In the LO state, there are two possible modulation direc- 
tions with respect to the applied magnetic field: longitudinal 
and transverse. In this letter, we consider the longitudinal LO 
state in the vortex lattice. Hereafter, the longitudinal LO vor- 
tex state will simply be called the FFLO state. To investigate 
the stability of the FFLO state against a magnetic field, we 
solve the microscopic Eilenberger equations self-consistently 
in the three-dimensional (3D) space composed of the in-plane 
vortex lattice and the longitudinal FFLO modulation, taking 
into account the orbital and Pauli-paramagnetic depairings on 



an equal footing. 

Furthermore, to provide fundamental theoretical informa- 
tion on physical quantities in the FFLO state, we examine the 
effects of the FFLO modulation on the nuclear magnetic res- 
onance (NMR) 9-12 ) and flux line lattice (FLL) form factors 
obtained in a small-angle neutron scattering (SANS) 13 ' 14) ex- 
periment. It will become apparent that these two methods can 
provide direct and crucial evidence of the FFLO state, among 
the variety of other experiments. 3, 15) We will discuss anoma- 
lous behaviors in the corresponding experimental data on the 
heavy Fermion superconductor CeCoIns, the high-field and 
low-temperature superconducting phase of which is regarded 
to be a realization of the FFLO state. Thus the main purpose 
of this letter is to demonstrate the NMR spectrum and FLL 
form factors through the theoretical study. 

Our basic strategy is to provide if -dependent properties of 
the FFLO states for the 3D Fermi sphere and s-wave pair- 
ing. The corresponding 3D calculation for the FF state 16) and 
the full self-consistent analytical theory for a quasi- ID case 8) 
have been performed previously in the Pauli limiting case 
without vortices. Here, we extend those calculations to take 
account of vortex effects. Before discussing the anomalous 
behavior of the FFLO states in CeCoIns, it is necessary to 
clarify the quantitative properties of the FFLO state in a typi- 
cal example of the 3D Fermi sphere. 

We calculate the spatial structure of the vortex lattice state 
using the quasi-classical Eilenberger theory in the clean limit 
valid for » 1 (It? is the Fermi wave number and £ is the 
superconducting coherence length). 17, 18) The Pauli paramag- 
netic effects are included through the Zeeman term jd B B(v), 
where B(r) is the flux density of an internal field and //b is 
a renormalized Bohr magneton. The quasi-classical Green's 
functions g(oj n + ifiB, k, r), f(oj n + ifiB, k, r), and f\co n + 
ijiB, k, r) are calculated in the vortex lattice state by the Eilen- 
berger equations 19, 20) 

{o^+ViB + v-(V + iA)}/ = Ag, 

{^+i/i5-v-(V-iA)}/ t =A*g, (1) 

where g = (1 - //^) 1/2 , Reg > 0, v = v/vfo, and the Pauli 
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parameter \i - \i^B^lnk^T z . k is the relative momentum of 
the Cooper pair, and r is the center-of-mass coordinate of the 
pair, v is the Fermi velocity and vfo = ( v2 )\J 2 where (■ • ■ >k 
indicates the Fermi surface average. We assume a magnetic 
field is applied to the z-axis. The Eilenberger units of Rq for 
lengths and Bq for a magnetic field are used. 19 ' 20) The order 
parameter A and the Matsubara frequency oj n are normalized 
in units of nk^T c . 

For self-consistent conditions, the order parameter is calcu- 
lated by 



A(r) = goN T Yj (f + f*\ 



(2) 



0<a) n <a) cl 



(3) 



with (g N ) 1 = InT + 2T Zo<^<^ cut a) n x . We use co cut = 
20&b7c- B = V x A is self-consistently determined by 

IT 

V x (V x A) = V x M para (r) - — ^ (vim g) k , 

K 0<cj„ 

where we consider both the diamagnetic contribution of su- 
percurrent in the last term and the contribution of the param- 
agnetic moment M para (r) = (0, 0, M para (r)) with 



M para (r) = M 



m 

H 



2T 
JlH 



0<co n 



(4) 



The normal state paramagnetic moment Mo = (jj/k) 2 H, k = 
BolnkQT c VSttTVo and No is the density of states at the Fermi 
energy in the normal state. We set the GL parameter k to be 
102. We solve eq. (1) and eqs. (2)-(4) alternately, and obtain 
self-consistent solutions, as in previous works, 19,20) under a 
given unit cell of the triangular vortex lattice. For the FFLO 
state, A(r) has a periodic oscillation with the period L in ad- 
dition to the vortex lattice structure. As the unit cell size of 
the vortex lattice is determined by H = (B) r , we can estimate 
the //-dependence of the FFLO state in our calculation of the 
vortex lattice. Throughout this paper we use /i = 5 as a repre- 
sentative case of the strong Pauli paramagnetic effect. 

The Gibbs free energies are calculated from self-consistent 
solutions using eq. (9) in ref. 21 for the FFLO state with var- 
ious FFLO wavelengths L. We compare them to identify the 
most stable state under a given H and T. Figure 1(a) exhibits 
the resulting successive changes at T/T c = 0.1. It is seen that, 
starting from H = H C 2 where the FFLO state with the shortest 
wavelength L = 17 is stabilized, L becomes longer as H de- 
creases. Eventually the free energy of the FFLO state becomes 
comparable to that of the conventional Abrikosov state, where 
the FFLO modulation along the field direction is absent. The 
envelope of the free energies of the FFLO state approaches 
that of the Abrikosov state, such that the two curves seem to 
merge tangentially, namely, at the meeting point, tangents of 
the two curves coincide with each other. While our calcula- 
tions are for discretized L, even these results suggest (1) a 
second-order-like transition between the FFLO state and the 
Abrikosov vortex state 7) and (2) the continuous L change of 
the FFLO state as a function of //, which are similar to the 
results of previous analytic FFLO theory. 8) 

In Fig. 1(b), we show the phase diagram in the H-T plane. 
This is obtained by repeating the FFLO calculations as a func- 
tion of H at different temperatures, T/T c =0.l, 0.15, 0.2, and 
0.25. H L q is the transition field from the Abrikosov vortex 
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Fig. 1. (Color online) (a) Free energy differences F from the normal state 
for the FFLO state with different wave numbers L and the Abrikosov state, 
as a function of H. T = 0.1T C . In the normal state, F = 0. (b) Phase 
diagram for the FFLO state in the H-T plane for the 3D Fermi sphere and 
the s-wave pairing. //= 5. H C 2 is the first order at H > H cr . Lines are guides 
for the eye. 
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Fig. 2. (Color online) Spatial variations of (a) the order parameter A(r) and 
(b) the paramagnetic moment M para (r) for several values of L along the 
field direction outside of vortices. These are normalized by its length and 
maximum values. T/T c = 0.1 and /u = 5. 



state to the FFLO state and the transition at H C 2 to the normal 
phase is the first order above H cr . The FFLO region in the H- 
T plane is given by //lo/#c2 = 0.973 for the present \i - 5 
at T/T c =0.l. The value of Huj/H C 2 depends on the \i value, 
namely, H L0 /H c2 = 0.991 for// = 2 at Z7Z>0.1. Even in this 
strong paramagnetic case of [i = 5, the FFLO phase appears 
only near H C 2, and //lo increases on lowering T in this typical 
example of an isotropic Fermi sphere. 

Figure 2 displays normalized waveforms of (a) the order 
parameter A(r) and (b) the paramagnetic moment M para (r) in 
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(a) 



Fig. 3. (Color online) Field evolutions of various quantities at T/T c = 0.1 
and jj, = 5. (a) FFLO wave number q = 2n/L. (b) Form factor |Fiool 2 - Inset 
shows the overall variation, (c) Form factor IF102I 2 . (d) Form factor IF104I 2 . 



the FFLO state along the field direction outside of vortices. 
It is seen that a simple sinusoidal modulation waveform con- 
tinuously deforms into an anti-phase kink form, or solitonic 
waveform as H approaches the #lo line where L diverges. 
In other words, near the //lo boundary, the sign change or 
71 phase shift of the order parameter occurs sharply, meaning 
that the excess electrons and M para (r) are confined in a narrow 
spatial region along the kink position. 

The FFLO nodal kink forms a sheet of paramagnetic mo- 
ments perpendicular to the field. On the other hand, along the 
vortex lines, enhanced M para (r) at the vortex core is decreased 
at the intersection with FFLO kink plane (see Fig. 2 in ref. 
19 and Fig.l in ref. 22). There, the zero-energy peak states 
of quasiparticles are absent, because of the 2n phase shift of 
the order parameter, coming from the kink and from the vor- 
tex. The paramagnetic moment becomes strongly confined to 
the kink position as H approaches //lo from above. As will 
be seen later, these three-dimensional FFLO spatial structures 
can be probed by SANS experiments or NMR experiments. 

As shown in Fig. 3(a), the FFLO wave number q = 2n/L 
of the stable FFLO state continuously varies with H . Starting 
with q = at H = H^o, q increases sharply. Hence the an- 
tiphase solitonic waveform quickly changes into a sinusoidal 
one on increasing H (see also Fig. 2). This behavior is similar 
to that seen in the exact solution (see Fig. 9 in ref. 8), implying 
that the FFLO physics along the parallel direction exemplified 
here is common and universal. 
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Fig. 4. (Color online) NMR spectra in the FFLO state: (a) paramagnetic 
moment distribution P(M) (b) and internal field distribution P(B). jd = 5 
and T/T c = 0.1. Upper panels show //-evolution of spectra in stereo- 
graphic view. Lower panels show spectra at some values of H. Horizontal 
baselines for each spectrum are shifted by H/H C 2, which is indicated on 
the right axis. 



The FLL form factor is an important quantity that can be 
directly measured in a SANS experiment. The form factors 
Fhu are Fourier components of an internal field B(r) in our 
calculation. 19) The fundamental Bragg spots F100 for the vor- 
tex lattice are shown in Fig. 3(b) as a function of H. The 
spot |Fioo I 2 increases in the Abrikosov state, because M para (r) 
accumulates at the vortex core to increase B(r) locally, as 
seen in the inset of Fig. 3(b). This feature has already been 
shown theoretically 20) and observed in various paramagnet- 
ically enhanced superconductors, such as TmNi2B2C 23) and 
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CeCoIri5. 14) As shown in Fig. 3(b), the intensity of l^iool 2 sud- 
denly decreases upon entering the FFLO phase and continues 
to drop quickly, almost exponentially. (Note the T = 50 mk 
data in Fig. 1 of ref. 24). This is because M pa ra(r) is not en- 
hanced at the vortex core on the FFLO nodal plane (see Fig. 
5(b) in ref. 19). This contribution decreases l^iool 2 , which is 
the average along the z-axis. 

In addition to the usual Bragg spots Fiqq associated with 
the vortex lattice, the observation of extra spots Fio n in =2,4, 
...) is crucial to prove the existence of the FFLO phase. In 
Fig. 3(c), we show IF102I 2 , which is the new superspot as- 
sociated with the FFLO modulation along the field direc- 
tion. It rises quickly at H = H^o- After attaining a maxi- 
mum in the middle of the FFLO phase, it slowly decreases 
towards H C 2. Thus the best chance to observe it is in the mid- 
dle field region inside the FFLO phase. The relative intensity 
l^i02l 2 /l^iool 2 = 1/10 ~ 1/20. Therefore it is quite possible 
to detect it because l^iool 2 is enhanced by the Pauli effect even 
near H c2 . The higher order spot I^kmI 2 is also shown in Fig. 
3(d). It takes a maximum just near H L0 . Since the magnitude 
of 1^104 1 2 is further reduced and is one order of magnitude 
smaller than IF102I 2 , it might be difficult to detect it. 

The NMR spectrum is also crucial to identify the FFLO 
state. By choosing probed nuclei that have different hyper- 
fine coupling constants, we can effectively pick up the se- 
lective field distributions. 19) When the hyperfine coupling is 
sufficiently strong, the paramagnetic distribution M para (r) is 
probed in NMR experiments. In contrast, in the weak cou- 
pling case, the magnetic induction B(r) in the whole system 
is detected by NMR. In the mixed state of ordinary supercon- 
ductors, it yields the so-called Redfield pattern. 

Here, we examine the field evolution of the NMR spectra 
for both strong and weak hyperfine coupling cases. For the 
former (latter), we evaluate the distribution P(M) [P(B)] using 
the stable FFLO state at each field. These are given by 

P(M) = J (M - M para (r)) dr, P(B) = J (B - B(r)) dr, (5) 

i.e., volume counting for each M and B. Figure 4 shows the 
spectral evolutions of these distributions for two cases. In Fig. 
4(a), P(M) is displayed. Since in the Abrikosov state the para- 
magnetic moment is confined exclusively to the vortex cores, 
a single main peak appears at the saddle point (S) position 
in the NMR spectrum. In the FFLO phase, M para (r), which 
comes from excess electrons, accumulates at the normal state 
(AT) position. The peak at the N-position becomes dominant 
towards H c2 , because an increasing excess of unpaired quasi- 
particles appear at the FFLO nodal sheets. It is noted that near 
H L0 , two peaks appear simultaneously in the NMR spectrum. 
This double peak structure was observed in the In(2) site of 
the NMR spectrum by Kumagai et al l2) for CeCoIns. The 
appearance of the double peaks at S- and N-positions is un- 
ambiguous evidence of the FFLO state. 

It is also important to observe the characteristic change of 
P(B) for the weak hyperfine case, exemplified by In(l) in 
CeCoIns. Here, rather unexpectedly, the double peak struc- 
ture can be seen in Fig. 4(b) near H^o, beyond which the N 
peak dominates the spectrum. The N-position is near the S- 
position in P(B), compared with P(M). In the lower field of 
the Abrikosov state, the usual Redfield pattern is reproduced, 



as seen from Fig. 4(b). 

We touch upon the recent NMR experiment on CeCoIn5. 12) 
The observed double peak structure of In(2a) for if || c and 
H || ab is markedly similar to our result in Fig. 4(a) (see the 
spectral evolution in Fig. 2 of ref. 12). The proposed phase 
diagram of FFLO for H \\ c is also similar to our Fig. 1(b) 
where #lo/#c2 ~ 0.975 for// = 5 compared with #lo/#c2 = 
4.77/4.95 T ~ 0.95 at zero temperature for H || c. lV} As men- 
tioned previously, the value of H^o /H C 2 depends on /d, but the 
shape of the FFLO phase diagram is hardly changed by the 
value of fi. For H || ab, the proposed phase diagram is modi- 
fied because of the presence of the existing SDW. 25) 

We also calculated the FFLO structure in d-wave pairing 
or the quasi-two-dimensional Fermi surface when a magnetic 
field is applied to the z-axis. These situations do not qualita- 
tively change the results shown in Figs. 2-4. Careful estima- 
tion of the free energy to determine the stable L remains for 
future work. 

In conclusion, we quantitatively explored the field evolu- 
tion of the FFLO state for typical examples of a Fermi sphere 
and s-wave pairing, by self-consistently solving the micro- 
scopic Eilenberger equations in the 3D space of the vortex 
lattice and the FFLO modulation along the field direction. To 
facilitate the identification of the FFLO state through experi- 
ments, we estimate the NMR spectrum and FLL form factors 
as a function of the magnetic field in the FFLO vortex states. 
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